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Abstract 

We discuss W and Z production through the deep inelastic viN- scattering 

in the context of the standard model SU(3)c x SU(2)l x U(l) of the strong 

and electroweak interactions. We find the cross section rates for the process 

vi+N — > l~ + W + + X for the case of ultrahigh-energy neutrinos (10 14 eV < 

E u < 10 17 eV) colliding on a target nucleon (Ej\f = mj\f). We also calculate 

a{y\ + M l~ + X) in order to compare it with a{yi + M — > l~ + W + + X). 
13.15.+g, 13.85.Tp, 25.30.Pt 
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I. INTRODUCTION 



Large-scale neutrino telescopes [1] have as a main goal the detection of ultrahigh-energy 
(UHE) cosmic neutrinos (E u > 10 12 eV) produced outside the atmosphere (neutrinos pro- 
duced by galactic cosmic rays interacting with interstellar gas, and extragalactic neutrinos) 
[2]. UHE neutrinos can be detected by observing long-range muons produced in charged- 
current neutrino-nucleon interactions. A very enlightening discussion on UHE neutrino 
interactions is given by R. Gandhi et al. [3]. 

The detection of UHE neutrinos will provide us with the possibility to observe vM- 
collisions with a neutrino energy in the range 10 12 eV < E u < 10 21 eV and a target nucleon 
at rest. Earlier estimations of the W production rates in z/A/"-scattering were given in 1970, 
by H. H. Chen [4], and by J. Reiff [5], in 1971 by R. W. Brown and J. Smith [6], and 
by F. A. Berends and G. B. West [7]. In those works the calculations were made using 
neither the standard model [8] nor the parton model [9]. One of the aims of detecting UHE 
neutrinos is to look for physics beyond [10] the standard model [8]. However, the event rates 
for most of the exotic processes could be of the same order as the event rates for heavy 
gauge boson production in z/A/"-scattering. Hence, even though the cross sections for heavy 
boson production are expected not to be large enough to allow for detailed investigations 
of the heavy bosons properties, it is necessary to calculate them. Only if heavy boson 
production, which could be an important background for physics beyond the standard model, 
is completely known and subtracted from the exotic event rates, the new physics could be 
investigated. Therefore, we discuss in this paper heavy vector boson production in the deep 
inelastic z^A^-scattering within the frame of the standard model SU(3)c x SU{2) L x U(l) 
of the strong and electroweak interactions and using the parton model [9] with the parton 
distribution functions reported by J. Pumplin et al. [11]. We use the CTEQ PDFs provided 
in a 71/ = 5 active flavours scheme. 

The cross section for the process v\ + M — > l~ + W + + X is expected to be the most 
important one because it gets contributions from photon exchange diagrams (see figs. 2.1- 
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2.5). We thus present results for the total cross section for W + production through the 
processes Vi+Af — > l~ + W + + X as a function of the neutrino energy in the range 10 14 eV < 
E v < 10 17 eV 1 and taking E_\f = m_\f for the nucleon. We also calculate a(vi+N — > l~ + X) 
to find out how it compares to a{yi + M — > Z~ + VF + + X). 

This paper is organized as follows. In section II, we provide the expressions for the 
kinematics of the production of a massive vector bosons in deep inelastic neutrino-nucleon 
scattering and discuss the influence of setting cuts on the phase space of the produced 
particles. In section III, we perform the explicit calculation of the matrix elements and the 
differential cross section for polarized scattering. We present and discuss our results for the 
total cross section for the case of W + production in section IV. Finally, section V contains 
a summary of our conclusions. 

II. KINEMATICS 

In this section we will discuss the kinematics and phase space of the production of a 
heavy vector boson B = (W ± , Z°) through the inclusive processes 

Vi+Af^l' + B + X (1) 

where u h J\f, I' and B stand for the incoming neutrino, the target nucleon, the final lepton 
[yi or l~) and the produced massive boson, respectively. We will denote the four-momenta 
of these particles by p, P_\f, p' and k, respectively. In accordance with the kinematics for the 
collision of a neutrino on a target nucleon, the following construction is chosen: 

p» = E u (l, 0, 0, 1), PjO- = T7w(l, 0, 0, 0). (2) 

As usual we define the invariants: 



x We take E v < 10 17 eV, in order to use the CTEQ PDFs in the range of the (a/, Q' 2 ) parameter 
space in which these have been found to be valid: 10~ 6 < x' < 1 and 1.3 GeV < \JQ' 2 < 10 4 GeV 
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and the dimensionless variables: 
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x 



2u r ' 



2v' 



(4) 



The physical region of these kinematical variables is obtained by requiring that the scalar 
products of any two particle four-momenta be positive, the determinants A 3 of any three par- 
ticle four-momenta positive, and the determinant A 4 of the four independent four-momenta 
(whenever possible we will neglect the nucleon and lepton masses): 



A 4 (P,P',P^ k) 



pp' pPj^f pk 

p'p p'Pm p'k 

Pmp Pnp' P N k 

kp kp' kP M M 2 B 



(5) 



negative [12]. 

From the positivity of the scalar products we find: 



< x < x' < 1, 



< y' < y < 1. 



(6) 



Explicit evaluation of A 4 using (3), (4) and (5) gives: 



where 



A 4 = -( S /2)VV (T 6 -r)(r-T a ) 



r aife = (1 - x')(l +y'-y) + (l/y){(l - y)[(x' - x)(y - y') - n] + xy' 
T2y/xy'(l - y)[{x> - x){y - y') -//]}. 



(7) 



(8) 



with /j, = Mg/s and s = 2m^E v . The condition A 4 < can be used to restrict r: 

T a <T< T b , (9) 

From the requirement that r a ,b be real, together with x > 0, y' > and y < 1, there follows 
the inequality 

(x'-aOfe-j/') >//, (10) 

which is equivalent to A 3 (p,p',k) > 0. Inequality (10) can be used to replace (6) by the 
more stringent relations: 

< x < 1 — u/y, x + iily < x' < 1, 

A* < y < 1, < y' < y-n/{x-x'). 

The expressions given in (8), (9) and (11) define the physical region for the dimensionless 
variables x, x 1 ', y, y' and r. We have taken all fermion masses to be zero, this implies that 
in the calculation of the total cross section the integration over the momentum transfer 
square extends up to zero. Hence, we have no natural regulator for the contribution of the 
photon-exchange diagrams. This is one of the reasons why we need cuts on Q 2 and Q' 2 . 
Another reason is that the parton distributions can be used only when Q 2 and Q' 2 are not 
too small. Furthermore, in order to separate deep inelastic from elastic scattering, a cut on 
the invariant mass W of the unobserved particles in the final state is required. Therefore, 
besides the kinematical conditions (8), (9), (11) we have in general also the following: 

Q 2 = sxy > Q 2 C , 

Q' 2 = sx'yi > Q'l (12) 
W = sy'(l-a/)+m%- > W c . 

The cuts for Q 2 , Q' 2 and W constrain further the physically allowed region for the process 
(1). This region can now be written in terms of the dimensionless variables as follows: 



with 

y a = max{(M B + sjw c ) 2 /s, (M 2 B + Q 2 c )/s}, y b = 1, 

x a = Q 2 Jsy, x b = l-(M B + ^W c ) 2 /sy, 

x' a = max{^[y(l +x) + (Mi - W c )/s - yj[y(l - x) - (M| + Wy/s] 2 - 4M' Wy S 2 ], 
^[xy + (Ml + g'2)/ a + ^ + (M 2 + g'2)/ a ]2 _ AxyQ'l/s}}, 

< = ^{^(1 + *) + (Ml - W c )/s + ~ x )~ i M l + W c )/s]> - AM^WJ? }, 

y' a = msx{W c /s(l - x'), Q ,2 Jsx'}, y' b = y- M 2 B /s(x' - x). 

The limits for r are not modified by the cuts given in (12), thus remaining as in (9). 

III. THE DIFFERENTIAL CROSS SECTION FOR INCLUSIVE BOSON 

PRODUCTION 

The differential cross section da UlAf for (1) is calculated in the parton model from the 
cross section da Viq of the parton subprocess 

vi + q -> I' + q' + B (14) 

(where vi, q, I', q' and B stand for the incoming neutrino, a quark inside the target nucleon, 
the final lepton {i>\ or the outgoing quark and the produced massive boson) and the 
parton distribution functions f q (x*, Q 2 ), which are the probabilities to find a quark q with the 
fraction x* of the nucleon momentum: g M = x*P t ^ 1 in a scattering process with momentum 
transfer square Q 2 . We will denote the four-momenta of these particles by p, q, p', q' and k, 
respectively. As usual, we define for the parton process the invariant variables s, Q 2 , ...,£>' 
and using the definitions given in (3) we get the following relations: 
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s = (p + q) 2 = x*s, 

Q 2 = -{p-p'f = Q\ 

v = q(p — p') = x*v, 

V ; (15) 

§' = (p + q-k) 2 = s' - (1- x*)s + 2(l- x*)(u -u'), 

C = q(p — p' — k) = x*v', 
The variables s, . . . , v' are not independent. For massless partons we have Q' 2 = 20' and 
consequently x* = Q' 2 /2v' '. Comparing with the definition of x' given in (4) we conclude 
that x' = x*, i. e. q» = x'Pfr. 

The parton cross section is obtained from the invariant matrix element M(uiq — > Vq'B): 

da viq = El-^l 2 dT 3- (16) 

The sum is performed over all fermion spin states s and the polarizations A of the produced 
boson. With help of the different sets of variables introduced in (15), (3) and (4), and by 
using (7), the 3-particle phase space dT 3 can be expressed as follows: 

d 3 p f d 3 q' d 3 k' 4 , , 

Zhipi Z,£liqi Z,£Lik 

7r dQ 2 d0ds'dQ' 2 d0' - /2 . 
= — — Q\Q —2v) 

» s J-A 4 (p,q,p',k) 



Ti dQ 2 du ds' dQ' 2 du' r/W2 „ , A 

o(Q — 2x u ) 



8s J-A 4 (p,P M ,p>,k) 
ns dx dy dy' dr 



(17) 



8 J(T b -T)(T-T a ) 

From (16) and (17) we obtain 

d^ii - 1 v i m^\ 2 - dxd y d y' dr 
vw# M 1 x' v / (T6 _ r)(T _ To) - 

The Feynman diagrams for .M^ 9 are depicted in fig. 1. The heavy boson can be produced 
from the lepton line (fig. 1: a, b), the (anti) quark line (fig. 1: c, d) or via the non-Abelian 
vertex from the line of the exchanged bosons (fig. 1: e). We write 

M Uiq = M Uiq ' 1 + M Uiq ' h + M Uiq ' n . 



Diagrams containing the exchange of Higgs bosons have been neglected because of the small- 
ness of the Higgs fermion coupling. We present the explicit expressions for these matrix 
elements in a way which is suited not only for the standard model of the electroweak inter- 
action but also for extended models containing more vector bosons than the SM. Therefore, 
denoting with / the fermions and with B the gauge bosons, we define the couplings: 

i> r tPfB^ i^(vf lf - af, /75 ) = i^(Lf lf P L + Rf, f P R ), 

and 

Kip) Kip') B p(p") i 9BB'B"{g^(p -p') P + 9u P (p' - p'% + g w (p" -p)u}, 
with Pl,r — (1 T 75)/2, L f, j = v f, f + af, f , and Rf,f = vf,f — af l f. 

In the case of the standard model the couplings of the fermions to the Z and W bosons are 
given as follows (e = \f^na): 

v f'f = ~ e ( T f ~ 2QfSin 2 6 w )/sin26 w , 
a z flf = -eTf/sin26 w , 
v J f = -e/(2V2~ sin9 w ), 
afif = —e/{2\/2sin9w)- 

Qf and T| stand for the charge and the third component of the isospin of the fermion f. For 
the three bosons non-Abelian couplings we have: 

9jW+w- — e > 
gz°w+w- — e cos6 w / sin6 w . 

The invariant matrix elements M UlQ,a , a = {I, h, n} for the production of the boson B (polar- 
ization vector E^(k, A)) in the scattering of the neutrino with initial and the outgoing lepton 
with final longitudinal polarization Pi, Pi and the quark with initial and final polarizations 
P q , P q i can be expressed, taking all fermion masses to be zero (hence Py = Pi and P q i = P q ) 
in the following form 

Mp g p = J h /« fp a P Tr j 1 / 1 
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where the left-, right-handed currents of the leptons and quarks are defined as usual 

4' fl =fi(p')7^«(p), 

and 

J^ p «=u(q') lfl P q u(p). 
and the tensors ', r = {1, 2, 3} are defined as follows: 

= = e"k u - k»e u , T% v = ie^ pa e p k a . (18) 

The invariant functions f P ^ P{ are built from the propagators and couplings in the form: 

/feU = 2e P Cplpt + 2e P' c p{,p fp^ = Cp^Pt ~ C p{,po 

fp l q , Pl = {-i)Hc l l Pl + c l l Pi \ 

Jp v ,p t ~ Z£( l p g ,p, + Z£ Q °p g ,Pi' Jp q ,Pi ~ u P q ,p t ~ u P q ,Pii 

fp q ,pi = 2£ (p-p') c p q \po fp q ,Pi = ~ c p q ,pv 

ft* = o, 

with Si = 0, 1, 0, 1 and S 2 = 0, 0, 1, 1 for P q Pi = LL, LR, RL, RR, respectively, and 
C l i] L = E Lf^L^/KQ" + Ml,){Ml - 2kp)} 

B',l" 

C l L f L = E LfrL^LfjKCr + Ml,){Ml + 2kp')} 

B',l" 

C^l = E LZ>Lf: q ,,L* q /[(Q* + M B ,)(M B - 2kq)} (19) 

B',q" 

C h L ' f L = E L%,L* q L*J[(Cf + ]tf B ,)(]tf B + 2k</)] 

B',q" 

C n L , L = E L^L^g B/BBII /[(Q 2 + M 2 B/ )(Q'" + M 2 Bf ,)}. 

B',B" 

for left-handed polarization of leptons and quarks. For other polarizations L has to be 
replaced by R in a suitable way in (19). 

The next step in the calculation of dap^ Pi is to square A4 P ^ Pi : 
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and sum over all fermion spin states s and the polarizations A of the produced boson. 
We have 



E^p^pS)* 



s,X 



Ejh,P q rr,a jl,P t ( jt,Pi\* fr ,b i ^ v \ */ A^y 

J /x iPq^P^r J v \ J v> ) JP^PjWr' J 1<V / 

s,\,r,r' 

E1_/Pq fPi £T,a <r-fivfr',b /,£•// 
ri f i,fi> L 'uv'JP q ,P l J ~r JP q ,P l \ J 'r' ) 

X,r,r' 



Erpr,r' rr,a fr',b 
I P q ,P l Jp q ,P l JP q ,P l 



with 



and 



J-h,P q n-h,P q 
II O u 



rn — <ri,Pi <ri,Pi 



rpr,r' _ n/Pq sppv ( spu'v' \* r>Pi 
1 P q ,Pl — rL niJ 1 ' J ~ r Wr> ) L -uv'- 



Using the expressions given for H^i and L p u i in the appendix and the fact that the / 



Pi 



are real functions, it is straightforward to show that 



M^(M^y = (M^yM^. 



Hence 



E \mZp\ 2 = E (\M^k 

s,X s,X 



+ 



MTp, 



+ 



M 



viq,n 
P q ,Pl 



+2Re ]T {M^iM^T + M^ l Pl (M^y + M^(M^y) 

{ s,X 

Explicit expressions for the quantities T r,r ' are presented in the appendix. Also, the explicit 
summation over the polarizations A of the produced boson is presented there. 

A heavy weak boson B can be produced in deep inelastic v\M collisions 2 via the following 
processes: 



2 the results for i^A^-scattering can be obtained from those of zV^-scattering by the replacements 
v\ — ► i>i, W ± — > W T , I' — > l + and u-type quarks <-> d-type quarks 
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vt + M + Z° + X, (P.l) 

v t +Af^l- + Z° + X, (P.2) 

v l +Af^v l + W+ + X, (P.3) 

W+N^>vi + W- + X, (P.4) 

u t +Af^l- + W+ + X, (P.5) 
ui+J^^l- + W~ + X, 



The diagrams which contribute in the lowest order at the quark level to the different reaction 
mechanisms (production at the leptonic vertex, at the hadronic vertex and through the 
boson self interaction) of all these processes, are depicted in figs. 2.1 - 2.5. The process 
vi + M — > l~ + W~ + X, is forbidden in the lowest order. We see from these figures that 
the reactions (P.l) - (P-4) only get contribution from heavy boson exchange diagrams and 
therefore their total cross sections are expected to be very small and we will not be discussed 
in this work. 

The final step in the evaluation of da UlAf consists now in putting together the parton cross 
sections da m and the parton distribution functions f q (x', Q 2 ). In contrast to deep inelastic 
iVV-scattering the choice of Q 2 is not unambiguous in our case since the momentum transfer 
to the nucleon depends on the reaction mechanism (in other words, whether the boson is 
emitted at the leptonic or at the hadronic vertex). For leptonic production it is reasonable to 
take Q 2 = — (p—p' — k) 2 = Q' 2 since p—p' — k is the momentum transfer to the nucleon. In the 
case of hadronic production the obvious choice is Q 2 = — {p—p') 2 = Q 2 ■ For the non-Abelian 
diagrams a simple kinematical argument is not sufficient. For heavy boson production in 
neutrino quark scattering, unitarity would be violated without a coupling between the W, 
Z and 7 bosons. Therefore unitarity is restored through strong cancellations between these 
non-Abelian diagrams and either the leptonic or the hadronic contributions to neutrino 
lepton scattering. 

In this work, we will neglect the contribution from heavy boson exchange diagrams. 
Hence, we calculate the cross section of process v\ + M — > l~ + W + + X with the following 
practical prescription which guarantees the restoration of unitarity 
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da^ = £ / dx'f q (x', Q' 2 ) (dar + daZ q + do™) . (20) 

q J 

This prescription collects the expressions where the heavy boson is emitted from the lepton 
line or the non-Abelian vertex and the interference of these reaction mechanisms. 

IV. RESULTS FOR W + PRODUCTION VIA vM COLLISIONS 

We have already pointed out in the previous section that the process v x + M — > l~ + W + + 
X is the only one which gets contribution from photon-exchange diagrams. Therefore, we 
will restrict ourselves to calculate numerically the cross section for H^ + -production through 
this process, leaving out the contribution from heavy boson exchange diagrams. 

We obtain our numerical results using the standard model of the electroweak interactions 
[8], taking M w = 80.4 GeV for the mass of the charged boson W and M z = 91.2 GeV for 
the mass of the neutral boson Z (hence, sin 2 9w = 0.223 for the electroweak mixing angle) 
[13]. We present results for the case of unpolarized deep inelastic z/A/"-scattering with a 
neutrino energy in the range 10 14 eV <E V < 10 14 eV and the nucleon at rest 1 , i.e. a target 
nucleon. We take cuts of 4 GeV 2 , 4 GeV 2 and 10 GeV 2 for Q 2 , Q' 2 and the invariant mass 
W, respectively. These values for the cuts are suited for the parton distribution functions 
of J. Pumplin et al. [11] which we will use in our calculations. In fig. 3, we show the 
dependence of the cross section for the reaction v x + P — > l~ + W + + X on the cuts on 
the momenta transfer square (Q 2 , Q' 2 ). Taking a cut W c = 10 GeV 2 for the invariant mass 
W, we plot a{yi + P -> l~ + W + + X) for the following cases: a) Q 2 C = Q'\ = 2 GeV 2 , b) 
Ql = Q'l = 4GeV 2 , c) Q\ = Q'l = 6 GeV 2 and d) Q 2 C = Q' 2 C = 8 GeV 2 . We observe in this 
graph that the total cross section rates of the process v\ +Af — > l~ + W + + X do not depend 
strongly on the choice of the cuts on the momenta transfer square (Q 2 and Q' 2 ), when they 
take on values of a few GeV 2 . 

In fig. 4 we plot o{y x + P — > l~ + W + + X) and a[yi + P — > l~ + X) as a function of 
E v . We obtain a{v x + P -> Z" + + X) = 2.3 x 10" 35 cm 2 for E v = 10 17 eV. In fig. 5, 
we show o[yi + N — > Z~ + + X) and cr(^ + A — > Z~ + X) as a function of £„. We get 
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a(vi+N -> 1-+W++X) = 1.8xl0~ 35 cm 2 for E v = 10 17 eV. Finally, in fig. 6, we present our 
results for a(u x +M -> /- + W^ + +X)/ ( r(z/ i +A/' -> Z~+X) (A/" : P, N) as a function of E v , with 
Pat = m/v. In particular, we obtain a(ui+P -> Z~+Py + +X)/o-(z/ ; +P -> Z~+X) = 4.7xl0 -3 
and cr(i/| + iV -> l~ + W + + X)/cr(^ + N -> Z~ + X) = 3.6 x 10~ 3 for P„ = 10 17 eV. 

V. CONCLUSIONS 

We have presented the general formulas for the cross section of the production of massive 
vector bosons in deep inelastic neutrino nucleon scattering in the framework of the standard 
model. The expressions for the matrix elements are given in a way that is suitable for 
extended models containing more vector bosons than the SM. 

We have neglected the contribution from heavy boson exchange diagrams in our numeri- 
cal calculations. Hence, we have given numerical results only for the total cross section rates 
of the process v x +Af — > l~ + W + + X (A/": P, N.) because this reaction is the only one which 
gets contribution from photon-exchange diagrams in the lowest order of a. Considering the 
nucleon at rest, taking M w = 80.4 GeV, sin 2 9 w = 0.223, a neutrino energy in the range 
10 14 eV < E u < 10 17 eV, setting cuts of 4 GeV 2 and 10 GeV 2 for the momenta transfer 
square (Q 2 and Q' 2 ) and the invariant mass (W). We made use of the parton distribution 
functions of J. Pumplin et al, we used the CTEQ PDFs provided in a nj = 5 active flavours 
scheme. In particular, we have obtained a{y\ + P — > l~ + W + + X) = 2.3 x 10~ 35 cm 2 and 
a{y x + N -> I- + W+ + X) = 1.8 x 10" 35 cm 2 for E v = 10 17 eV. We have also shown that 
the total cross section rates of the process v x + N — > l~ + W + + X do not depend strongly 
on the choice of the cuts on the momenta transfer square (Q 2 and Q' 2 ), when we take them 
equal to a few GeV 2 . 

We have also presented results for o{y x + M -> l~ + W + + X)/a(u x + N -> /~ + X) 
(Af : P, N) as a function of E v in the range 10 14 e\/ < E v < 10 17 e\/, with E M = m M . We 
have gotten a(v x + P ^ l~ + W + + X)/a(u x + P -> Z" + X) = 4.7 x 10~ 3 and (t(^ + X -> 
Z- + + X)/(j(^ + X -> Z- + X) = 3.6 x 10- 3 for E v = 10 17 eV. 
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APPENDIX A 



By setting all fermion masses to zero, in this appendix we give the expressions for the 
quantities T r ' r ', which are defined as follows: 

rpr,r' _ n/Pq <rix v ( T^' v ' \* T P l 
1 P q ,P l — rL /i f i'' r r Wr> ) L -uu'i 

with 



(r,r' = 1,2,3), 



Klfi = 2{q^> + ~ 9w>qq' T ie^ lpu q p q' a } 

j=r = 

T pv = £ »k v - e v W 

J 3 — it tph, a 

and P q , Pi being the polarization of the initial quark and initial neutrino, respectively. In 
general, these quantities are functions of scalar products of the momenta p, q, p', q' , k and 
the polarization vector e of the produced boson. 

Using the polarization sum for the massive vector boson 



J2 A)e„(fc, A) = -g^ + 



(Al) 



and the definitions given in (18) we get 

Y, JTi&TY = -9^ + 9^ k v k p + 9 UP kilk ° - 9™ kPkP 

A 
A 

£ JTW)* = ^ P 9 va - 9 iUJ 9 vp )M 2 B + E-FfW)* 

A A 

Aided by these expressions we obtain 

8[-M|(pp' • qq' + p'q' ■ pq - p'q ■ q'p) 
+2{pq • kp' ■ kq' + p'q' • kp • kq)} 



A 



22 



LL, RR 



8[—Ml(pp' ■ qq' — p'q' ■ pq + p'q ■ q'p) 
+2(pq' ■ kp' ■ kq + p'q ■ kp ■ kq')] 



LR, RL 
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Y^ReT 

A 



2.3 



J2 T 

A 



33 



±16 (—pq ■ kp' ■ kq' + p'q' ■ kp • kq) 

=Fl6(p'g • kp ■ kq' — pq' ■ kp' ■ kq) 

8[Ml(pp' • qq' — p'q ■ pq' — pq ■ p'q') 
+2(pq ■ kp' ■ kq' + p'q' ■ kp ■ kq)] 

8[Mg(pp' ■ qq' — p'q ■ pq' — pq ■ p'q') 
+2(pq' ■ kp' ■ kq + p'q ■ kp ■ kq')] 
The remaining T,s can be expressed as follows 

IQp'q' . pq 

16p'q ■ pq' 

—8(ep ■ p'q' ■ qk + ep' ■ pq ■ q'k 
—eq ■ p'q' • pk — eq' • pq ■ p'k) 



T ll = { 



ReT 



12 



ReT 



13 



— 8(ep ■ p'q ■ q'k + ep' ■ pq' ■ qk 
—eq ■ pq' ■ p'k — eq' ■ p'q ■ pk 

=F8(ep • p'q' ■ qk — ep' ■ pq ■ q'k 
+eq' ■ pq ■ p'k — eq ■ p'q' ■ pk) 

±8(ep ■ p'q ■ q'k — ep' ■ pq' ■ qk 
—eq' ■ p'q ■ pk + eq ■ pq' ■ p'k) 



Using (Al) and (18) we get 



Y j eP l -eP 2 = -P l P 2 + 



kP 1 ■ kP 2 



Ml 



Y^eP-n v = H V {e^-P) 

A 

5>P • = J=r(e -> -P) 

A 



for P, Pi, P 2 = p, q, p', q'. Hence 

J2ePi-eP 2 T 

A 



11 



-PlP2 



kP 1 -kP 2 . 

Ml ■ 



LL, RR 

LR, RL 
LL, RR 

LR, RL 



LL : RR 

LR, RL 
LL, RR 

LR, RL 
LL, RR 

LR, RL 



T 



11 



16 



Y,eP-ReT 12 = ReT l2 { 

A 

J2^P-ReT 13 = ReT l3 ( 



17 



Figure Captions 



Fig. 1: Feynman diagrams for heavy boson production from the initial (a) and final (b) 
lepton, from the initial (c) and final (d) quark and via the non-Abelian couplings (e). 
Figs. 2.1 - 2.5: Feynman diagrams for processes (P.l) - (P-5): boson production from 
the incoming neutrino (a), the outgoing lepton (b), the initial (c) and final (d) quark, and 
through the non-Abelian couplings (e). u stands for u, c, d, s, b; d for d, s, b, u, c. 

Fig. 3: a{y x + P -> l~ + W + + X) as a function of E u in the range 10 14 eV < E u < 10 17 eV, 
with Ep = m P . Taking W c = 10 GeV 2 and a) Q 2 C = Q'\ = 2 GeV 2 (curve A), b) Q\ = Q'\ = 
AGeV 2 (curve B), c) Q 2 C = Q' 2 C = 6 GeV 2 (curve C), d) Q 2 C = Q'\ = 8 GeV 2 (curve D). 

Fig. 4: a[yi + P — > l~ + W + + X) and a(^ + P — >■ /~ + X) as a function of in the range 
10 14 e\/ < ^ < 10 17 e\/, with E P = m P . Taking W c = 10 GeV 2 and = Q'2 = 4Ge y2 > 

Fig. 5: o[yi + X — > l~ + M 7 " 1 " + X) and cr(i/j + X — > l~ + X) as a function of E v in the range 
10 14 eV <E V < 10 17 e y, with E N = m N . Taking W c = 10 GeV 2 and Q\ = Q'\ = 4 GeV 2 . 

Fig. 6: cr(^ + M -> Z" + + X)/a(^ + AT -> Z" + X) (A" : P, N) as a function of 
E v in the range 10 14 eV < E v < 10 l7 eV, with E N = m N . Taking W c = 10 GeV 2 and 
Q 2 c = Q' 2 c =AGeV 2 . 
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